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ABSTRACT 


To assist analysis of riveted lap joints, stress intensity factors are determined for surface and 
comer cracks emanating from a wedge-loaded hole by using a 3-D weight function method in 
conjunction with a 3-D finite element method. A stress intensity factor equation for surface cracks is 
also developed to provide a closed-form solution. The equation covers commonly-encountered 
geometrical ranges and retains high accuracy over the entire range. 


INTRODUCTION 


Current damage tolerance analysis of aircraft structures is mainly based on linear elastic fracture 
mechanics, where the stress intensity factor is a key parameter. Surface and comer cracks emanating 
from circular holes are among the most common defects occurring in aircraft structures. Stress 
intensity factor solutions for these cracks are available in the literature for remote tension and remote 
bending [1-5]. Based on three-dimensional finite element solutions, the corresponding stress intensity 
factor equations were also developed [3], and are being used in a variety of applications. However, 
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for another important category of loading conditions (i.e. wedge loading of a cracked hole) very 
limited solutions [1,4] exist in the literature for comer cracks, and no solution is available for surface 
cracks. The primary interest in the wedge loading case lies in the fact that solutions for commonly 
encountered pin-loading case, such as the riveted lap joints shown in Fig.l, can be obtained by 
superposing solutions for remote tension, remote bending and wedge loading. The objective of this 
paper is to provide stress intensity factor solutions for surface and comer cracks emanating from a 
wedge-loaded hole, which can be used to obtain solutions for rivet loading. This is accomplished in 
two ways. First, an accurate and efficient 3-D weight function method [5] is extended to solve the 
problems by using stress distributions from a 3-D finite element analysis of the same, but otherwise 
uncracked, configuration. The weight-function analysis code covers a wide range of geometrical 
parameters and various loading conditions, and produces accurate solutions on Personal Computers. 
Second, a stress intensity factor equation for surface cracks is developed by curve-fitting procedures. 
This latter effort is to facilitate the determination of stress intensity factors for commonly-encountered 
surface cracks in damage tolerance analysis. 


METHODS OF ANALYSIS 


The crack configurations to be considered are shown in Fig. 2 together with the definition of 
relevant parameters (infinite plate width and height assumed). The 3-D weight function method is 
used in combination with the 3-D finite element method to solve these problems. This combined 
approach takes optimal advantage of both methods. The weight function method is accurate and 
efficient in solving 3-D crack configurations. However, its accuracy can not be guaranteed without 
using, as input, an accurate stress solution for the same, but otherwise uncracked, configuration. For 
crack configurations involving stress concentrations (such as a hole or notch), the stress distribution in 
the region of uncracked configuration where cracks will occur has a fully 3-D nature, the stress varies 
not only in plate- width direction but also in thickness directions, even though the applied load does 
not change through the plate thickness. This 3-D effect manifests itself most near the intersection of 
the hole with the plate surface for small r/t ratios (see for example [6,7]). Therefore, to account for 
this 3-D nature, a 3-D finite element analysis was performed [7] to provide the normal stress 
distributions. Because the finite element analysis involves no cracks, only one finite element 
calculation is required for each r/t ratio. The normal stress distribution obtained from the finite 
element analysis is then fitted into the following equation: 

a(x,y)/a 0 = EEcyipW* (1 > 

«=l j - 1 r 1 

Equation (1) is used in the following weight function equations: 
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K<M>a x ) =j° x [o(x,;y)-Pc*,;y)] W a (a x ,y)dy 


(2a) 


K c<y>Cy> - f 0 Cy P(x,y)W c (c r x)dx 


(2b) 


where K.fraJ and K c (y, Cy ) are, respectively, stress intensity factors for a-slice and c-slice, intersecting 
at point fry) on the crack front (see Fig. 2 (d) and (e)), and W a and W c are the weight functions. The 
stress intensity factors for the 3-D crack, K(<p), are then obtained by eq.(3)- 


K(ip) = 


(- 1 )" 


l-ri 2 (v,a/c,Aq>) 


K*(x,a x ) 


Y K e(y>c y ) 


(3) 


For further details of the weight function method, please refer to [5]. 


RESULTS AND COMPARISONS 


The stress intensity factors are presented in a dimensionless form defined as 


F(<p) =K(<p)/(oy/na/Q) 


(4) 


in which a is a reference stress, and equal to remote uniform stress for remote tension, and a=P/(2rT) 
for wedge loading, where P is the total applied force along the bore of the hole. 


Corner Cracks 


Reference [1] provides stress intensity factors for wedge loading with a cosine squared pressure 
istribution, while the present weight function solutions are for wedge loading with a cosine pressure 
distribution. However, the difference caused by the different pressure distributions will be negligible 
as long as the c/r ratio of cracks exceeds about 0.8. The following comparisons are made for such 
c*ses. Figure 3 shows the comparison of the weight function results with the finite element solution 
[1J- Very good agreement is observed. Figure 4 (a) shows the difference between a single comer crack 
and double comer cracks. Figure 4 (b) gives the stress intensity factor ratio of double comer 
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cracks over single comer cracks. The corresponding estimate by using Shah’s empirical formula [8] 
and the ratios for through-the-thickness cracks [9] are also shown in Fig. 4 (b). 


Surface Cracks 


The comparison with finite element solution [3] for surface cracks is made for remote tension, 
since no solution is available for surface cracks under wedge loading in the literature. Figure 5 shows 
the comparison. Again, very good agreement is observed. To examine the behavior of the weight 
function method, some limiting cases are also considered as a/c-^O and oo. Figures 6 (a) and (b) show 
these limiting cases. It is clear that the appropriate 2-D limits are reached under the limiting 
conditions. 


STRESS INTENSITY FACTOR EQUATION FOR SURFACE CRACKS 


In many situations, a closed-form solution is desirable, as evidenced by the popularity of the 
previous stress intensity factor equations for other cases [3]. Therefore, an effort is made to develop 
an equation for double surface cracks under wedge loading based on the weight function solutions. 
The applicable range of the equation is chosen to be 1 <r/t<5, 0.5 <a/c<2 and 0.005 <a/t< 0.9. 
T his range is believed to cover most cases encountered in practice. If a problem does go beyond this 
range, then the weight function code is directly used. The equation is developed separately at 0 and 
90° only, so as to reduce the number of independent variables from 4 to 3 (r/t, a/c and a/t). The 
general form of the equation is: 

K(<p) = ay/na/QF(<p) <p =0° ,90° (5 

The developments for F(90°) and F(0°) are described in the following sections. 


Development of F(90°) 


The form of F(90°) is: 
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(6) 


F(90°) =(-) r f 9 (-,— ,—) 
c t c t 


where 7 — 0 for a/c < 1 , and 7 =-! for a/c> 1. The f 9 function is given below: 


f/ r a a* f r a a., ,x a. 

/ 9 ( 7 ,-, 7 )-*, ( 7 ,-, 7 ) /„( 7 , 7 ) 


(7) 


where g 9 is a fitting function, and f sa is a combination of the two functions f s (x/r) and f a (a/t). 

The f s (x/r) is the dimensionless normal stress distribution along the x-axis for a wedge loaded hole in 
an infihite plate. According to Rooke & Tweed [10], f s (x/r) is expressed as: 


/&. 2 S &.1 

ran 


3 +- 


1 + 2(1 +— ) 2 - 3 (l +— ) 4 
r r 


( 1 +-) 

r 


hJ 


2(1 +-) 3 

r 


2 +— 


x 

r 


( 8 ) 


The f a (a/t) is the dimensionless stress intensity factor for collinear cracks, and given as [11]: 



tanC-K-) 


ll 


(9) 


The introduction of f s and f a is based on the consideration that F(90°)=f s (0)*f a (a/t) as a/c^0. Using f s 
and f a , f sa is defined as: 




3/,(0)/.(5) +U-VJ0) 


(10) 


The use of f sa is to bring in some of the contributions of the rest of the crack to F(90°). The use of 
these functions is believed to be very helpful in fitting complicated functions involving multiple 
variables in a wide range. The g 9 is fitted in the following form: 




(ii) 


The functions a 1 (r/t,a/t) and a 2 (r/t,a/t) are expressed as: 
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“•<? f ' =f ■* ( 7> +#1 4 ) * ( 7 )3 } * exp b < 7 ) * ( ?> 


,/*„ ,r “m 

a ii(-^) = a m *(-) 


« isC^:> = « 131 + « 132 * c^ } + a 133 ( Y } 


“ m( 7*7 


,r a. f ,/\ .a. ,7\ .a. 2 l 

“ 2 ( 7 . 7 ) ={“2,(7) + “22( 7 )*( 7 )*“23( 7 )*( 7 ) }*< 

“2.< 7 )"“2U + “2i2*( 7 ) + “20*( 7 ) 2 + “2«*( 7 ) 3 . i-1.2.3.4- 


The a parameters with 3-digit subscripts are constants, and listed in Table 1. 


( 12 ) 


(13) 


Development of F(0°) 


The functional form of F(0°) is given as: 


F(0°) = 

c r c r 


(14) 


where X=l/2 for a/csl, and X=- 1/2 for a/c>l. 'The f 0 function is given below: 

/o Vc? ' go Y?7 )/c V 0 Y?7 ; 


(15) 


where g 0 and t 0 are fitting functions, and f c is a weighted average of f^c/r). The latter 

is the dimensionless stress intensity factor for through-the-thickness cracks under the same loading, 

and given as: 

fj£.) = 1.788 -4.935(-) +8.527(-) 2 -7.790(-) 3 +3.450(-) 4 -0.5832(-) 5 (16) 

cu r r r r r r 

Then, f c is defined as: 
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(17) 



Similarly to F(90°), the introduction of is based on the consideration that F(0°)=f c0 (c/r) as 

a/c->oo . And f c is used to bring in some of the contributions of the rest of the crack to F(0°). The go 
is fitted in the following form: ^ 


t c t t t c 


(18) 


The functions ftCr/t.a/t) and 0 2 (r/t,a/t) are expressed as: 



Pi 2 (f)*(f) + Pi 3 (f)*(f) 3 


P 1 1(~) “ P 1 11 + P 1 12 *(-) 

Pl2(“) = Pl21 + Pl22*(T) + Pl23*(~) + Pl24*(“) 
* t t t 

P 13(~) = P 131 + P 132 * (“) + P 133(~) 


(19) 



The (3 parameters with 3-digit subscripts are constants and listed in Table 2. The ^ is a fine-tuning 
function and equal to 1 for a/t<0.2, and fitted as follows for a/t>0.2: 
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(21) 


^ ,r a a s ,r a. . ,r a s ,<*,.4, ,r a Xj ./a. 2 

o( 7 ~c =jPi(t ’ t> jP2( "’ “ } ( “ } Ps( “’ T> ( “ } 


r a, 
> > 
f t 


t t 


t t c 

3 


p,(-.-)=p»(-) + p u (-)*(T) + Pa(r)*<-:) i=1 - 2 - 3 - 


r t 


pA i- o. 12,3 

l Z t 


The constants p ijk for to are given in Table 3. 


Comparison of Equation with Weight Function Data 


The accuracy of the equation is examined by comparing it with the original weight function data. 
Figures 7 (a) and (b) give F(90°) and F(0°) versus a/t, respectively. The variations of F(90°) and 
F(0°) with a/c are shown in Figs. 8 (a) and (b), respectively, which are for a/t=0.9, the largest a/t 
ratio considered. The largest deviation of the equation with the weight function data occurs at this a/t 
ratio for F(90°). From Fig. 8 (a) the maximum deviation is within 3%. Figures 9 (a) and (b) show 
F(90°) and F(0°) versus r/t, respectively. It is seen that the equation is generally within 2% of the 
weight function data, with maximum error of about 3% in the whole range considered. 


CONCLUSIONS 


Based on the above analysis and results, the following conclusions can be made: 

1 . The combination of the 3-D weight function method and the 3-D finite element method is an 
efficient approach to solve 3-D crack problems involving stress concentrations. 

2. Accurate stress intensity factors are provided for surface and comer cracks emanating from a 
wedge-loaded hole. 

3. A stress intensity factor equation is developed for double surface cracks. The equation covers 
commonly-encountered geometrical ranges and retains very good accuracy. 
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Table 1. Coefficients of g 9 function 



1 

2 

3 

4 

am 

«l2i 

a i3i 

1.292 

-0.6158 

-0.3834 

-0.02280 

0.03094 

0.1790 

-0.02065 


a2ii 

a 22i 

a 23i 

a24i 

0.06285 

0.1827 

0.09046 

1.871 

0.01986 

-0.2219 

-0.02106 

-0.7457 

-0.008287 

0.08763 

-0.02337 

0.4400 

0.0009832 

-0.01045 

0.005125 

-0.05090 


Table 2. Coefficients of go function 


1 

2 

3 

4 

Pin 

0.9708 

-0.005293 



&\2i 

-1.415 

0.7752 

-0.1815 

0.01502 

013i 

0.4926 

-0.1570 

0.01652 


0 21i 

0.0007527 

0.05418 

-0.02149 

0.002344 

&22\ 

0.4667 

-0.4266 

0.1586 

-0.01709 

&23i 

3.345 

0.04976 

-0.6747 

0.1018 

&24\ 

-3.626 

1.941 

-0.2248 




Table 3. 

Coefficients of ^ function 


1 

2 

3 

PlOi 

6.89367 

-53.6151 

50.8248 

Plli 

-4.08538 

40.9391 

-46.3484 

Pl2i 

0.677014 

-9.35989 

12.1718 

Pl3i 

-0.0226694 

0.668142 

-1.00152 

P20i 

-12.0523 

108.472 

-102.842 

p21i 

8.63812 

-84.3581 

94.9457 

p22i 

-1.60061 

20.0463 

-25.5194 

p23i 

0.0778192 

-1.51170 

2.15760 

p30i 

5.83664 

-53.1766 

50.5411 

p31i 

-4.20676 

41.6555 

-46.9805 

p32i 

0.801768 

-10.0724 

12.7830 

P33i 

-0.0421292 

0.778241 

-1.09613 
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Fig. 1 Riveted lap joints with cracks. Fig.2 (a) wedge-loaded hole, (b) surfece cracks, (c) corner cracks, 

(d) & (e) crack front coordinates. 



K d /K 



Fig. 3 Comparison of dimensionless stress intensity 
factors with finite element solutions. 


Fig. 4 (a) Comparison between single crack and 
double cracks. 




Fig. 4 (b) Ratios of stress intensity factors of 
double cracks to single cracks. 


Fig. 5 Comparison with finite element solutions for 
surface cracks under remote tension. 
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Fig. 6 Limiting behavior as (a) a/c tends to 0, 
(b) a/c tends to infinity. 
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Double surface cracks 
Wedge loadinq (cos) 
a/t=0.9 


1.5 


Double surface cracks 
Wedge loading (cos) 
a/t = 0.9 




Fig. 8 Comparison of the equation with original weight function data, (a) at 90 degree, (b) at 0 degree. 




Fig. 9 Dimensionless stress intensity factor 


as a function of r/t ratio, (a) at 90 degree, (b) at 0 degree. 



